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ABSTRACT 



Context. The first observational evidence of a connection between supemovae (SNe) and gamma-ray bursts (GRBs) was found about 
a decade ago. Since then, only half a dozen spectroscopically confirmed associations have been discovered and XRF 1003 16D asso- 
ciated with the type-Ic SN 2010bh is among the latest. 

Aims. We constrain the progenitor radius, the host-galaxy extinction, and the physical parameters of the explosion of XRF 
1003 16D/SN 2010bh at z = 0.059. We study the SN brightness and colours in the context of GRB-SNe. 

Methods. We began observations with the Gamma-Ray burst Optical and Near-infrared Detector (GROND) 12 hours after the GRB 
trigger and continued until 80 days after the burst. GROND provided excellent photometric data in six filter bands covering a wave- 
length range from approximately 350 to 1800 nm, significantly expanding the pre-existing data set for this event. Combining GROND 
and Swift data, the early broad-band spectral energy distribution (SED) is modelled with a blackbody and afterglow component atten- 
uated by dust and gas absorption. The temperature and radius evolution of the thermal component are analysed and combined with 
earlier measurements available from the literature. Templates of SN 1998bw are fitted to the SN itself to directly compare the light- 
curve properties. Finally, a two-component parametrised model is fitted to the quasi-bolometric light curve, which delivers physical 
parameters of the explosion. 

Results. The best-fit models to the broad-band SEDs imply moderate reddening along the line of sight through the host galaxy 
(Av.hoit = 1.2 ± 0.1 mag). Furthermore, the parameters of the blackbody component reveal a cooling envelope at an apparent initial 
radius of 7 x 10" cm, which is compatible with a dense wind surrounding a Wolf-Rayet star. A multicolour comparison shows that SN 
2010bh is 60-70% as bright as SN 1998bw. Reaching maximum brightness at 8-9 days after the burst in the blue bands, SN 2010bh 
proves to be the most rapidly evolving GRB-SNe to date. Modelling of the quasi-bolometric light curve yields Mn] = 0.21 ± 0.03Mo 
and Mj.j = 2.6 ± 0.2Mq, typical of values within the GRB-SN population. The kinetic energy is E^^ = (2.4 ± 0.7) x 10^^ erg, which is 
making this SN the second most energetic GRB-SN after SN 1998bw. 

Conclusions. This supernova has one of the earliest peaks ever recorded and thereafter fades more rapidly than other GRB-SNe, 
hypemovae, or typical type-Ic SNe. This could be explained by a thin envelope expanding at very high velocities, which is therefore 
unable to retain the y-rays that would prolong the duration of the SN event. 

Key words, gamma-ray burst:individual:XRF 100316D - supemovae:individual:SN 2010bh 



^ 1. Introduction 



The core collapse of massive stars is thought to give rise to su- 
pernovae (SNe) and lon g-duration gamma-ray bu rsts (GRBs, du- 



rations larger than 2 s: iKouveliotou et alj|1993 ). The first clue 



to the connection betwe en these events is the sim ilarity in their 
kinetic-energy scale (see Wooslev & Bloomll2006l for a review). 
Even before the discoveries of the late 90's, a few authors dis- 
cussed the possible association between high-energy outbursts 
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and SNe dColgatJ 1 968t IPaczvriskill 1 986l) . In theory, the collaps- 
ing core of a very massive star, whose envelope has been blown 
away by its own stellar winds, i.e., a Wolf-Rayet (WR) star, 
can lead to the formati on of a relativi s tic jet that will produce 
high- energy emission ( Wooslev 1993 : Wooslev & MacFadvenI 
19991) in the form of a GRB or an X-ray flash (XRF softer 
events th ought to be mainly produced by the same mechanism 
as GRBs: lHeise et alJbOOlUKippen et al.ll2004l: Isakamoto et af 



2005':'Heise & in 't Zandl200ll) . WheiTthe jet colHdes with the 
circumstellar mate rial, it produces a multi-w avelength afterglow 
(for a review, see IZhang & MeszarosI 120041) . Theoretically, the 
energy of the core collapse should also be capable of expand- 
ing the remaining envelope. However, it is unclear how this can 
happen, or even if there is always enough energy for a SN ex- 
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plosion. The kind of SN thought to be associated with GRBs are 
those called stripped-envelope (SE) SNe, whose hydrogen enve- 
lope have mostly been removed. Events such as these are unusu- 
ally energetic stellar explosi ons and sometimes labelled "hy per- 
novae" in the literature (e.g., |Paczvnskilll998HHansenlll999l) . 

To date, only type-Ic SNe have been related to GRBs. The 
first and most representative case was that of SN 1998bw, 
whic h was found to be associated with the so ft GRB 980425 
(e.g., iGalama et all Il998bt iKippen et all [l998i) . A couple of 



days a f ter the trigger SN 1 998bw was discovered (iGalama et al 



1998at ISadler et al. 1998) inside the 8' error circle of GRB 



980425 (iSoffitta et all 11998.) in the underluminou s late-type 
galaxy ES0184-G82 (z = 0.0085 : | Tinnev et a J 19981). Al though 
initially controversial ( Galama etalj 1998bt ' 



Plan et al. 199 



the physical association between these objects was supported on 
temporal and spatial ground s by the slowly variable X-ray source 



at the position of the SN ( Plan et al. 2000l: Kouvehotou et al 



20041) . As a result of this connection, both the GRB and SN re- 
search fields, which until then had both evolved more or less 
independently, were revolutionised by a single event. 

Five years after SN 1998bw, the association between GRB 



scopic identification (IHiorth et alJl2003; 


Kawabata et al.l 


2003; 


Stanek et al. 


2003 




Matheson et all 12003 


) and became the first 



truly solid piece of evidence in favour of the GRB-SN connec- 
tion. Moreover, after 1998 there have been many other spec- 
troscopic associations, such as the cases of G RBs 02121 1 (SN 



20021t; Delia Valle et al. 2003), 020903 (Soderbe rg etalJl2005 . 
Bersieret alfcood) . 031 203 (SN 20031w; |Malesani etaklbool) . 



050525A (SN 2005nc; iDella VaUe et all 2006b), 060218 (SN 



2006a i; .Pian etal. 20061) . 81007 (SN 2008hw: Delia Valle et al 
2008h. 091127(SN 2009nz: lBerger et al.l201 lUCobb et al.l20l'or 



Filgasetal.1 boilh . and 101219B (SN 2010ma: I Sparreetal 
201 1) based on hypernova features in their spectra 



Furthermore, late-time bumps in the light curves of GRB 
afterglow s have b e en iri te rpreted as SN signa l s, e.g., GRBs 
970228 dReichartl 1 19991: IGalama etal J l2000t Ifieichart et al 



2000l) . 980326 dCastro-Tirado & Gorosabellll999[ 

!< MM 11 P - 



Bloom et al 

T 



1999h. 0111 21 (Bloom et al. 2002; Greiner et al. 2003), 020405 
dPrice et al.E003: .Masetti et al...2003.) . 040 924 (.Soderberg et al 



200( 



Wiersema et al.l 



l2008h . 



and 090618 (both in lCano et al. 



041006 (Stanek etal 



]2QQ± 



Soderberg et aljT2006h. 050824 (ISoUerman et alJ^OO Y). 060729. 



201 lal) to mention a few. These 



bumps show consistency in terms of colour, timing, and bright- 
ness with those expected for the GRB-SN population, but they 
are usually faint, which hampers the spectroscopic identification. 
These re-brightenings ha ve been detected in GRB light curves 
out to redshifts of ~ 1 ('Masetti et al. 2005; ' Delia Valle et al 



2003t iBloom et al.ll2009: Tanvir et al.,,201 0) owing to the sensi- 
tivity of current ground-based telescopes dedicated to follow-up 
observations. Sample studies of GRB-SNe (including bumps not 
spectroscopically identified) have been carried out to determine 
the luminosity distribution, the morphology of the light curves, 
and the physical parameters of the explosiori such as kinetic en- 
ergy, ejected m ass, and ^ ^Ni mass ( Zeh et al]|2004 : Ferrero et alj 
2006HRichard son 2009i: lThone et alj|201 lb . It has been asserted 
that GRB-SNe are in general brighter than t he local sample of SE 
SNe, except for cases such as SN 2010ay (ISanders et alJl201 ll) . 



In contrast to the claims of IStanek et alJ (120051) . no correlation 
has been found between the brightness at maximum with the 
shape of the hght curve tFerrero et alji2006) . 

Additional information about the explosion can be obtained 
from its early emission. One particular case is that of the 



soft XRF 06021 8 asso c iated with SN 2006ai (ICamparia et al 
2006'; 'Pian et al.' '2006': iFerrero et al. 2006; Cobb et al. 200^ 



iModiaz et al...2006: .Sollerman et al..,2006: .Mirabal et al.„200i 
This displayed an early X-ray and UV emission, which 
was interpreted as thermal radiation produced by the shock 



breakout f rom the surface of the progenitor ('Colgatel 
Fau3ll978t ISlein & Chevalier 1978; Matzner & McKee 



1974 



1999; 



Waxman et alj 120071: iNakar & Saril 12010) . The envelope is 
heated up and owing to expansion it starts to adiabatically cool, 
which then shifts the emission to UV and optical wavelengths. 
From the analysis of this a signal, it is possible to constrain 
the afterglow component, derive both the temperature and lu- 
minosity of the thermal co mponent, and comp ute the appar- 



ent radius of emission (e.g., iThone et al.M2011l) . Other exam- 



npu 

a. 



pies include SNe without detectable y-ray emission that never- 
theless exhibit adiabatic c ooling in the UV/optical and/or X-ray 



Malesani et al 



obser v ations: SN 2008D dSoderberg et al.l[2008l: 
'2009': 'Modia z et all 12009^, 'Mazzali et al.' '2008*), SN 2008ax 
(Romingetal.l l2009l). SN LS-04D2dc ( Schawinski et al. 20081), 
and SN 2010aq dOeziriet al.ll2010h . to mention a few. As addi- 
tional information, it has been claimed that high-energy emis- 
sion in GRB-SNe comes from acceler ated shock-breako ut pho- 
tons rather than highly relativistic jets dWang et alJi2007l) . 



Whilst no SN signature is expected for short GRBs, which 
are thought to be produced by the mergers of compact objects 
(for deep no n-detectio ns, see the c ases of GRB 050509B in 
Hiorth et alj|2005a andiBloom et al.ll2006 . and GRB 050709 in 



Fox et al.l l2005l and iHiorth et al.l l2005a) . there are a few sup 



posedly long events where an expected SN appearance was 



Ofek et al.ll2007l) and 060614 ( 


Fvnbo et al.ll2006; Gal- Yam et al. 


2006; Delia Valle et al. 2006a), there are verv tight constraints 



on the SN signature, which go down to 1% as bright as SN 
1998bw. The validity of the non-detections of SNe components 
for these two GR Bs as a classification tool is a point of con- 
troversy (e.g., Zhang et alj|2009) . We refer to Kann et al. ( 2011 ) 



for an exhaustive discussion on GRBs 060505 and 060614, in 
the light of the optical luminosity of their afterglows. 

We analysed the optical and near-infrared (NIR) data of XRF 
100316D and its associated SN 2010bh. The paper is organised 
by summarising the observations, data acquisition, reduction, 
and analysis in Sect.|2l The main results are treated separately as 
three different sections. The modelling of the early broad-band 
spectral energy distribution (SED) provides the progenitor radius 
and the host-galaxy extinction and is presented in Sect. [3] The 
multi-wavelength light and colour curves are analysed in Sect. |4] 
along with comparisons with previous GRB-SN events. In Sect. 
|5] the quasi-bolometric light curve is analysed and the physical 
parameters of the explosion are derived. Finally we gather our 
conclusions in Sect.|6l 
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,HST F814W, to+8.7d 




GROND i' image subtracted 
to+ 10.5 d 



N 



Fig. 1. HST/WFC3 F814W, GROND / '-band, and GROND host-subtracted / '-band images of tlie GRB-SN field, respectively, from left to right. 
Images were taken at around maximum brightness. Each panel is approximately 13" x 13" in size. The HST image shows significant galaxy 
structure near the explosion site (marked with two lines), which is blended with the object of interest in the GROND image shown in the middle 
panel. GROND images in the middle and right panels are shown using the same flux scale. 



2. Observations 

Gamma-ray emission fr om XRf 100316D trig gered the Burst 
Alert Telesco pe (BAT; Barthelrnv et al.l 120051) on board the 



Swift satellite (Gehrels et al. 2004) on March 16, 2010, at fo = 
12:44:50 UT (Stamatikos e t al.n201Q.). It t urned out to show a 
soft y-ray spectrum (Sakam oto et a n [2010h and a duration of at 



least 1300 s, one of the longest ever measured (iFan et al.ll201 1 



Starling et al.ll201 11) . About 15 h thereafter, a spec troscopic red 



shift o f 0.059 was published for the host galaxy (IVergani et al 



2010 a.b). Observations of the Gam ma-Ray burs t Opti c al and 



Near-infrared Detector (GROND; iGreiner et al.1 l2007l 120081) 



confirmed th at the new sotirce b ecame evident about 16 h af- 
ter the burst ( Afonso et al. 2010h . The rising of t h e sup ernova 



was verified photometrically by IWiersema et al.l (12010) only 



three days after the trig ger The spectroscopic confirmation of 
Chornock et al.l (1201 Obi) came approximately six days after the 



burst, w hich was confirmed two days later by iBufano et al 



(I2010ah . The S N was officially named SN 2010dh eight day s af- 
ter the trigger jBufano et al.ll2010bl: Ehornock et aklboiO c). On 
March 26, 2010, additional GROND observations of SN 2010bh 
were reported alon g with results a t the first attempts of host- 
galaxy subtraction (IRau et al ] l2010h . 



2.1. GROND 



The multi-channel imager GROND dOreiner et aDl2007l Eoosl) . 
mounted at the MPG/ESO 2.2m telescope on La Silla, Chile, 
started observations of XRF 100316D 11.7 h after the trigger, 
simultaneously in g'r'i 'z 'JHK^, with an average seeing of 1 '.' 1 , 
as soon as the astronomical night began. 

The GROND data were reduced using standard 
pyraf/IRAlf] tasks (^Todv 19931). s imilar to the procedure 



described in detail by Krti hler et al.l (l2008h . A general point- 



' IRAF, the Image Reduction and Analysis Facility, is distributed by 
the National Optical Astronomy Observatory, which is operated by the 
Association of Universities for Research in Astronomy (AURA), Inc., 
under cooperative agreement with the National Science Foundation 
(NSF); see http://iraf.noao.edUj 



spread function (PSF) model was constructed using bright field 
stars, from which the full width at half maximum (FWHM) was 
derived for each image. Aperture photometry was performed 
for science and calibration objects using an aperture size equal 
to the FWHM. The SDSS field at coordinates R.A.(J2000) 
= 06''59'"33^6, Dec.(J2000) = -17°27'00" was observed 
during photometric conditions to calibrate our images by 
performing relative photometry. A total of six stars in the field 
of XRF 1003 16D were employed for this purpose (see Table 
lA.ll of the online material). The same set of stars was used 
to calibrate JHK^ against the 2MASS catalogue. Calibration 
uncertainties vary in the range 0.002-0.020 mag for g'r'i'z' 
and 0.02-0.12 mag for JHK^, which together with catalogue 
systematics are added in quadrature to the statistical error 

A deep host-galaxy observation was carried out on 
November 5, 2010. This observation resulted in images with 
mean seeing of 0'.'74 and was used as a reference image for 
the subtraction of the host contribution from the early epochs. 
To develop a notion of the host contamination, we refer to 
Fig. \T\ where the position of the transient relative to the host 
galaxy is shown along with an F814W image of the Wield 
Field Camera 3 (WFC3) on board the Hubble Space Telescope 
(HST) as reference. In addition, the coordinates of the cen- 
tre of the host galaxy, catalogue name Anon J071031-5615, 
were obtained from the last GROND observation: R.A.(J2000) 
= 7^'10'"30^37(+0^07), Dec.(J2000) = -56°15'20'.'2(+0'.'3). 
To align the input and reference images, we use the WCSREMAP 
packag^ For the main purpose of subtracting template from sci- 
ence images, the HOTPANTS packag^l was employed. It matches 
the point spread function (PSF) and count flux of both input im- 
ages. It uses Gaussian functions to model the PSF in sub-regions 
of the original image. Since reduced images are essentially the 
combination of at least four g 'r 'i 'z ' and 24 JHKfi dithered in- 
dividual exposures, three and five Gaussian functions of differ- 
ent widths and degrees of freedom are employed for g'r'i'z' 
and JHK^, respectively. Point sources define the transformation 



- http://www.astro.washington.edu/users/becker/wcsremap.html 
^ |http://www .as tro.washington.edu/users/becker/hotpants.html 
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that is used to scale the template image to the flux of the sci- 
ence image. The routine outputs a noise map of the resulting 
difference image, which is employed to derive the uncertainties 
in the measured fluxes. The package DAOPHOT/IRAF was used 
to compute the magnitudes and their corresponding errors. This 
procedure was executed for a total of 140 individual images in a 
total of 20 epochs (see Fig. [1] for an example of image subtrac- 
tion in the / 'band). The position of the transient is R.A.(J2000) 
= 7^'10"30^55(+0^05), Dec.(J2000) = -56°15'20'.'0(+0'.'2) in 
host-subtracted optical images. The resulting photometry is tab- 
ulated in TablelA.2lof the online material. 



2.2. Sw\H/XRTand UVOT 

On board the Swjf t satellite, the X-Ray Telesc ope (XRT; 
Burrows et al.l2005b and the UVOT dRoming et al.l2 005) started 



observations of XRF 100316D at to + 2.4 min ("Stamati kos et aP 
[20 IQ) . Whilst a bright X-ray source was detected inside the BAT 
error circle, initial ly no afterglow candidate was found by UVOT 
dOates et al.ll2010l) . 

However, in deeper images taken at fo + 33 ks in the uvwl 
filter and at fo + 63 ks in the u band, we found evidence of emis- 
sion in excess of the host-galaxy contribution. To remove the 
contribution from the host galaxy, we requested ToO observa- 
tions in the uvwl and u filters, which were taken at fo + 3 x 10^ s 
(347 d after the burst), and amounted to a total exposure time of 
1525 and 1369 s, respectively. We measured a host galaxy con- 
tribution within the source aperture of 25 + 2 and 39 + 3 juJy 
in the uvwl and u bands, respectively. Subtracting this contribu- 
tion from our earlier-time data gave us a 3cr detection in 1894 s 
of uvwl- and in 4901 s of M-band data taken at mid-times of 
fo + 33 ks and fo + 63 ks respectively. In our analysis, we in- 
clude only the uvwl detection at fo + 33 ks, since there were no 
GROND data contemporaneous with the epoch of our M-band 
detection. A ll UVOT data were reduced following the procedure 
described inl Poole et al. ( 2008 ). To minimise the contamination 
from the underlying host galaxy, the source flux was measured 
within a circular source-extraction region of a 3 '.'5 radius. An 
aperture correction was then applied in order to remain compat- 
ible with the UVOT effectiv e area calibration s, which are based 
on 5" aperture photometry dPoole et al.ll2008r) . The background 
was taken from a source-free region with a 15" radius close to 
our source. 

The relatively bright X-ray afterglow (30-40 cnt s"' be- 
tween 144 and 737 s after the trigger) faded cons iderably at the 



begin ning of the second XRT epoch (fo -1- 33 ks; see lStarling et al 



201 li for a detailed analysis). In the subsequent analysis, we em- 



ployed XRT data at stages contemporaneous to GROND obser- 
vations, specifically in the interval from 33 to 508 ks after the 
burst. These data were obtained from the public Swift archive 
and reduced in the standard manner using the xrtpipeline 
task from the HEAsoft package, with response matrices from 
the most recent CALDB release. All data were obtained in photon 
counting mo de and downloade d from the XRT light curve repos- 
itory (E vans et alj|2007 , l2009h . Spectra were grouped using the 
grppha task. 

All the data discussed throughout the paper were corrected 
for the Galactic-foreground extinction of E(B - y)Gai = 0.117 



mag with Ry - 3.08 (ISchlegel et alJll998l) . All uncertainties in 
the following analysis are quoted at the Icr confidence level. 



3. Early broad-band SED 

Using data between 33 and 54 ks (roughly from 1 1 to 15 h) af- 
ter the burst, we compiled two early broad-band SEDs of XRF 
1003 16D. GROND provides detections in g'r'i'z'J and up- 
per limits in HK^, whilst the count rate of contemporaneous 
SwiftlXKY observations is already 0.01 cnt s"' at fo + 33 ks and 
decaying. After combining data from 33 to 508 ks after the trig- 
ger, XRT provides only three bins in the X-ray energy range. As 
follows, XRT data are scaled to the GROND first two epochs by 
using a decay index of a - -1.3 + 0.2, which is derived from 
the same XRT data over time. ^wZ/f/UVOT observed only in the 
uvwl filter at these stages. 

3.1. Modelling scheme 

Early blue emission coming from the XRF position was detected 
by GROND. For data acquired 42.5 ks after the burst, we mea- 
sured a colour g' - r' - -0.30 + 0.06 mag, in contrast to the 
red afterglows that usually follow a GRB (photon in dex F in the 
range 1 .2-2.5). Similar observations were made bv ICano et al. 
(1201 Ibl) . who found that this emission is incompatible with syn- 
chrotron radiation. The adiabatic cooling of the expanding atmo- 
sphere following the shock breakout gives us a reasonable ex pla- 
nation of the observed blue colours (e.g., ICano etal.ll20ilbh .In 
this scenario, the emission from the shock breakout lasts only 
a few hours after the core collapse and its SED resembles a 
blackbody at a high temperature of the order of 10^ K (or 0.1 
keV equivalently; e.g., Campan a et al. 2006h . Thus, we inter- 



pret the blue colours in our observations as the thermal compo- 
nent associated with the cooling of the shock breakout, which is 
similar to what has been claimed for XRF 060218/SN 2006aj 
( Campana et al. 2006; Waxman et al. 20 07|) and other ear ly- 
caught SN e (e.g.Tls oderberg et al. 2008; Modiaz et aDl2009l on 
SN 2008D; lRomin 'g et al. 2009 on SN 2008ax). In contrast to the 
idea of a thermal component producing the observed emission, 
there are no significant contemporaneous detections of UVOT 
in the uvwl bandpass, which alludes to high reddening. On the 
other hand, i'-J - 0.06 + 0.15 mag measured at the same epoch 
appears to have an underlying additive red component, which is 
represented in this case by the GRB afterglow synchrotron emis- 
sion (see Fig. IS . 

To test this hypothesis, we modelled the SEDs at 42.5 and 
50.0 ks (the latter shown in Fig. |2]i using two additive compo- 
nents: (1) a power law for the afterglow, and (2) an ideal black- 
body for the thermal component, respectively, of the forms 



PE(r)^CiE- 
Be(Tbb) 



and 

CiE'^dE 



(kTsB )4(ei^(i+:)/*7-BB - 1)' 



(1) 
(2) 



where Ci and C2 are normalisations, E is the spectral energy 
in units of keV, F is the photon index, k is the Boltzmann con- 
stant, Tbb th e intrinsic blackbody te mperature, and z - 0.059 is 
the redshift ( Chornock et al. 2010c). The model also accounts 
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Table 1. Fits to broad-band spectra constructed using GROND and SwiftjXKY data. 



Time Interval 
[s] after the trigger 


Reddening 
Law 


V)host 
[mag] 


kTYiB 
[eV] 


^BB 

[lO^^ergs^i] 


^BB 

[10'2 cm] 


r 


A^H.host 

[10^2 cm-2] 




42182-42879 
42182-42879 
42182-42879 
42182-42879 


MW 
LMC 
SMC 
MW 


0.21+oi| 

-0.21 

16+°" 

'-'■^"-0.16 

0.39 fixed 


-4 

81+*; 

82!^ 
79 + 3 


4.1+2-0 

-2.5 

4 O+I? 
3.4+1.7 
7.4 + 0.9 


27+8 
26+8 
24+^ 
39+4 


1.82+0.05 
1.82+0.05 
1.82+0.05 
1.73 + 0.06 


4.2+0.5 
4.2+0.5 

4.1+0.5 

4 4+0.5 


4.9/4 
4 9/4 

-r. yj-r 

ASIA 
8.7/5 


46630-53807 


MW 


0.39 ± 0.03 


78 + 2 


7.2 + 1.1 


38 + 3 


1.77 + 0.05 


4.4 + 0.4 


13/5 


46630-53807 


LMC 


0.38 + 0.03 


78 + 2 


7.0+1.0 


38 + 3 


1.76 + 0.06 


4.4 + 0.4 


16/5 


46630-53807 


SMC 


0.37 + 0.04 


79 + 2 


6.2 + 0.9 


35 + 3 


1.77 + 0.05 


4.4 + 0.4 


18/5 



Notes. The model consists of a blackbody plus a power law, both attenuated by optical/NIR reddening and X-ray metal absorption. 




100316D is employed as computed by the NED database^ fol- 
lo wing the standard ACD M model, using the redshift measured 



0.1 

Energy [keV] 

Fig. 2. Broad-band SED at 50 ks after trigger. The observed data are 
represented by black filled circles. The HK^ bands provide only 3cr 
upper limits shown as arrows. The thick grey line shows the best-fit 
model able to reproduce the data: extinguished power law plus black- 
body components, which are shown individually with dotted and dashed 
grey lines, respectively. The thin continuous grey lines show the unex- 
tinguished versions of the power law corresponding to the afterglow and 
the thermal component. In the lower panel, the residuals of the best fit 
are plotted. 



for host-galaxy extinction based on either MiUcy Way (MW; 
Rv = 3.08), Large Magellanic Cloud (LMC; Ry = 3.16), or 
Small Magellanic Cloud (SMC; Ry - 2.98) extinction laws and 
soft X-ray metal absorption. The Galacti c metal absorption i s 
fixed to be A^H.Cai = 7.05 x lO^o cm'^ dKalberla et al.ll2005h . 
The model has a total of six free parameters: C\, C2, F, Tbb, 
E{B - y)host, and A^H.host- The luminosity of the blackbody is 
computed as Lbb = 8.0525 €2(1 + z)D^^^ erg s"', where Dio is 
the luminosity distance to the transient in units of 10 kpc. A lu- 
minosity distance of 240 + 17 Mpc to the host galaxy of XRF 



by IChornock et al. ( 20I()cl) . a Hubble constant of 74.2 km s ' 
Mpc"' (lRiessetal.112009'), and the model of the local velocity 
field by Mould et al. (2000). 

The different set of parameters are summarised in Table [T] 
All best-fit parameters are consistent between the two epochs 
within their statistical uncertainties between both epochs. The 
temperature T^b and luminosity Lbb of the thermal component 
are in the range 78-81 eV and 4-7 x lO'*^ erg s"' from Cols. 4 
and 5 of Table [1] respectively. The afterglow power-law pho- 
ton index F a; 1.8 is shown in Col. 7. The reduced ;i^2 (or 
Xfj^ = X^/fi, where /i is the number of degrees of freedom) im- 
proves from 16.6 for the model without extinction to 2.7 for the 
model extinguished by MW-like dust. The best fit to the sec- 
ond SED epoch is shown in Fig. |2]with a thick grey line. The 
UV dust feature characteristic of the MW extinction law gives 
the more precise results, although it is poorly constrained at the 
bluer end. The X-ray tail of the blackbody fits the two data bins 
at around 1 keV, whilst the power law fits the only data point 
at X! 3 keV; both components contain significant absorption. We 
note that the fitting of the X-ray data has practically no residuals, 
i.e., the model over-predicts the data in this energy range. 

In addition to the luminosity and temperature, it is possible to 
compute an apparent emission radius of the thermal component 
(Rbb) from these two measured quantities by assuming isotropic 
radiation from an ideal blackbody. For all our trial models, the 

corresponding radii were calculated as /?bb = (4^'^^bb ^bb-* 
where cr is the Stefan-Boltzmann constant (Col. 6 of Table [T]i. 
Our best-fit parameters yield a radius of 3-4 x 10'^ cm, which 
is two orders of magnitude larger than the typical sizes o f the 
most likely GRB progenitors (WR stars; ICappa et al ] |2004 . 

The amount of metal absorption (estimated based on an 
equivalent hydrogen column density at solar metallicity) and 
dust extinction required for a good fit are A^H.host ~ 4 x 10^^ 
cm"2 and E{B - V )host - 0.2 - 0.4 mag, respectively. These two 
parameters are fitted independently and with no assumption be- 
ing made about the environment gas-to-dust ratio. However, the 
available data is not enough to allow us to distinguish among 
the different extinction laws employed in the modelling proce- 
dure. Whilst the reddening is consistent with values found in 

* |http://nedwww.ipac.caltech.edu/| 
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previous studies ( Starling et alJ2011 : Cano et al.l2011b ). the hy- 
drogen column dens ity differs significantly from the results of 



Starling et alj (1201 11) . The discrepancy is due to the additional 



constraints provided by the optical/NIR data, which were in- 
cluded to derive the parameters of the blackbody contribution. 
With our data set it is possible to tie the thermal component 
at both low and high energies, which provides a more accurate 
value for the absorption by heavy elements because of the greater 
constraint on the X-ray flux. 

3.2. Host-galaxy extinction 

We then attempted to evaluate the constraint on the extinction 
value derived in the previous section. Contour plots of redden- 
ing against blackbody temperature, luminosity, and spectral in- 
dex from the best fit to the data are shown in Fig.|3] respectively. 
The Icr and 2cr contours show that the solution to our best-fit 
relation is well-defined in all three panels. In the first panel there 
is a second solution to our best-fit relation at the 3cr level, which 
is at lower temperatures and similar reddening, although its sig- 
nificance is rather low. In the middle panel, the contours show 
a slight trend of proportionality between luminosity and extinc- 
tion, which is expected. At luminosities higher than 1.1 x 10"*** 
erg s"^ there is no possible solution, because the model becomes 
much brighter than the X-ray data and A^n.host cannot be allowed 
to vary and compensate for the X-ray luminosity of the model. 
Contours of the variations in the reddening and spectral index are 
shown in the lower panel, where the extinction is again restricted 
to high values. 

Analyses of optical spectroscopy of the host galaxy have 
found that the line-flux ratio of Ha to Hfi in the Hii region 
coinc ident with the SN position is consistent w ith zero extinc- 



2011; Levesque et al. 201 l l). Nevertheless, 



tion (IStarhng et al. 
th e Ay.host ~ estimate fr om the spatially resolved spectroscopy 
of iLevesque et al.l (1201 ih probes a much larger region 1.3 



kpc ) than the one probed by our line of sight to the GRB-SNe. 
A sufficiently high dust dumpiness could explain our high ex- 
tinction values along the line of sight and A i/ host ~ when in- 
tegrated over a larger patch. The position at which the hydrogen 
lines are formed, i.e., the H n region, might be located in front of 
the explosion site and not probe the same line of sight through 
the host galaxy. 

Alternatively a broken power-law model was fitted to the 
data in Fig. |2] for which Xt^^ - 3.0, which is slightly larger 
than that of the model that consists of a blackbody plus power 
law. The broken power law provided a closer fit when we as- 
sume that there is no host-galaxy extinction, significantly differ- 
ent from our results for the blackbody plus power-law model. 
However, the low-energy spectral slope of (3 - H-0.5 is incom- 



atible with synchrotron radiation (ySmax = +1/3; e.g.. lSari et al 



pane 



1998 ). Furthermore, the spectral break lies between the g ' and 
the uvwl bands at Vbreak = (8 + 3) x 10^ GHz, which is incon- 
sistent with the self-absorptio n feature usually observed at radio 
frequencies (vg ~ 2 - 13 GHz:lGalama et al. 19 98ctlTavlor et al 



1998t iGranot et al.1 119991; iGalamaet al...2000i) . Given also that 



Av,host > 0.2 mag is derived in Sect. 13. H and preferred by other 
authors (see next paragraph), the broken power-law model can 
be discarded with confidence. 



> 
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Fig. 3. Contour plots for the best-fit parameters from the modelling of 
the second epoch. From the outer- to the innermost contours, the green, 
red, and black lines are 3, 2, and Icr contours, respectively. The tail at 
low r in the bottom panel is caused by numerical inaccuracies. 



Additional evidence of large redden ing along the line o f 
sight inside the host galaxy was found bv lStarling et al.l (1201 ll) . 
E{B - y)host ~ 0.9 mag. The reddening was derived by fit- 
ting ^wi/f/BAT-i-XRT data and a M-band 3cr upper limit pro- 
vided by 5w!/f/UVOT data in the interval from 638 to 737 s after 
the burst. We employed a model that consists of a blackbody 
plus power law exting uished by SMC dust. A similar method 
was used here and by ICampana et al.l (120061) . who determined 
E{B - y)host = 0.20 for SN 2006aj. Another a ttempt at dered- 
dening SN 2010bh was carried out bv lCano et al.. (.201 lb) , who 
found E{B - y)host - 0.18 + 0.08 mag by assuming that the 
colours of SB SNe ar e all the same ten days after the V-band 
maximum brightness ( Drout et al. 2010h . Whilst this method is 
supported for the hydrogen atmos pheres of type-IIP SNe by 
a line of physical arguments (e.g., ^ Olivares E. et al. 2010l) . it 
is entirely empirical for SE SNe. Nevertheless, their redden- 
ing value is larger than zero with a significance of 2.3cr and 
consistent with our calculations for our second-epoch SED of 
E(B - y )host = 0.39 + 0.03 mag at the 2.5cr confidence level. 

After correcting for E(B - V )host = 0.39 + 0.03 mag, we 
obtained uvwl = 18.15 + 0.27 and u = 19.70 + 0.45 mag in 
the AB system at 33 and 63 ks after the burst, respectively. For 
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SN 2006aj ( Campana et al. 20061) . these values were uvwl - 
17.73 ± 0.21 and u - 17.77 + 0.15 mag at the same phase and 
redshift of SN 2010bh. If we had assumed that we should have 
seen a cooling envelope for SN 2010bh of comparable bright- 
ness and evolution as that shown by SN 2006aj, the host-galaxy 
extinction should have been higher than that estimated by our 
method. Nevertheless, the uncertainties are large and the com- 
parison of the uvwl measurements is consistent at the 1 .2cr con- 
fidence level. Although in the u band the uncertainty is even 
larger, there is no consistency with the brightness of SN 2006aj 
at the 4cr confidence level and a higher host-galaxy extinction 
('^y^host ~ 2.1 mag) would be required to reach the Icr level of 
consistency. 

In conclusion, we use the extinction value from the fit to 
the second-epoch broad-band SED throughout the paper, E(B - 
V\ost: = 0.39 + 0.03 mag for MW-like dust with Ry = 3.08, 
given that compared to the first epoch the statistical errors in the 
g 'r 'i 'z'J photometry are smaller and additional uvwl photom- 
etry is available. Moreover, when fixing the host-galaxy redden- 
ing to E(B - y)host = 0.39 mag, the quality of the fit to the 
first epoch is still acceptable (Xm^ - l-^)^ whilst the fit to the 
second epoch using E(B - y)host - 0.21 mag from the first- 
epoch modelling results vn - 3.7 and unphysical parame- 
ters. Hence, values of host-galaxy extinction including the cor- 
rection for redshift (^T-correction based on the spectral model) 
for the GROND filters and their corresponding statistical un- 
certainty are A^-.host = 1-50 + 0.12, Ar'.host = 1-10 + 0.09, 
A,',host = 0.80 + 0.08, Ae-host = 0.60 + 0.06, Ay,host = 0.39 + 0.04, 
Aff,host = 0.22 + 0.02, and A/f^.host = 0.14 + 0.01, all in units of 
magnitude. 

3.3. Evolution of the thermal component 

Having determined temperature, luminosity, and radius for the 
thermal component, it was then possible to study their evolution. 
Starhng et al.l (1201 ih analysed combined BAT+XRT data until 
737 s after the trigger and derived blackbody temperatures, data 
that we used in the following analysis. Temperature and radius 
over time are shown in Fig. |4] The decrease in temperature and 
the increase in radius are both trends that are consistent with 
the cooling of the envelope due to expansion and, thus, with the 
explosive scenario. 

The temperature and radius of an expanding envelope that is 
cooUng adi abatically have been theoretically shown t o evolve as 
power laws ( Wax man et al.l2007l:lNakar & Saril2 010V We there- 
fore fit power laws to the decay and rise of temp erature and ra 
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Fig. 4. Temperature and radius evolution of the blackbody compo- 
nent. Black dat a points a re those derived using combined BAT+XRT 
data (Starling et al.ll20llh plus the temperature determination s using 
GROND data (» 46 ks). The first data bin from IStarling et all j201lh 
represents the interval from -175 to 144 s, however, due to fitting and 
plotting purposes, here it corresponds to 0-144 s. Solid, dashed and 
dotted lines are different power-law models (see main text). Best fits 
to temperature and radius are shown in grey and brown solid lines, re- 
spectively. Grey crosses (upper panel) and open circles ( lower panel)) 
are measurements of XRF 060218/SN 2006aj taken from lKaneko et al] 
( l2007i) and iCampana et al.. (i2006i) . respectively. 



ofRo = (7.0 + 0.9) X 10" cm with Xti^ = 1-1 (solid line in 
the lower panel of Fig. |4]l. The resulti ng radius is slightl y larger 
than the size of WR stars (~ lO" cm: 'Cappa et al.l l2004V which 
are thought to be the pr ogenitors of long-du ration gamma-ray 
bursts and type-Ic SNe (IWooslev et al.ll2002h . Because of this, 
the initial emission radius of the thermal component might indi- 
cate the position at which a preexisting de nse wind surrounding 
the progenitor becomes opticallv thin (e.g..lCampana et al ■120061: 



diusto the combined data presented here and in IStarling et al 
(1201 iV The expression (0 = Ti — Kt was fitted to the evolv- 
ing temperature (solid line in the upper panel of Fig. |4]l, which 
gave a decay index of S - 0.3 + 0.2 for an initial temperature of 
r, = 0.17 + 0.04 keV withx/i ^ - 0.8, where k is the normalisa- 
tion of the power law. 

A model of the form ^BB(f ) = + vt^ was employed to 
fit the radius measurements. In the case of linear expansion, i.e., 
7=1, the radius grows at a mean velocity of v ~ 8,000 km 
s ' between BAT+XRT (until 737 s) and GROND observations the blackbody luminosity changes to 4-6 x 10*^ erg s"' in our 



ISoderberg et al.ll2008l:lBalbe rg & Loeb'201 iV From the theoret 
ical point of view, it is y = 0.8 (W axman et al...2007) . although, 
the data do not favour this solution delivering Xm^ ~ ^ (dotted 
line in the lower panel of Fig. |4]l. 

For an assumption of adiabatic cooling, the luminosity must 
be constant and was fixed to Lbb = 3.5 x l O'*^ erg s~', value de - 
rived from early-time X-ray measurements (IStarling et al .1201 ih . 
Assuming the best-fit model for the radius evolution from the 
previous paragraph, we computed the evolution of the adiabatic 
temperature as Tbb = (4;rcrLgg 7?^^ )"'/"', which is shown by the 
dashed line in the upper panel of Fig. |4] This model is clearly 
inconsistent with our data set; however, it indeed shows consis- 
tency with the XRF 060218/SN 20 06aj data set (grey crosses 
in the upper panel of Fig. |4l from iKaneko et al.l l2007b . Since 



(at 42-54 ks after the burst). When y was allowed to vary, we 
obtained a growth index of y = 1.4 + 0.3 and an initial radius 



late-time measurements, the assumption of a constant luminos- 
ity is invalid and inconsistent with the late-time temperature de- 
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termination. This result implies that either the cooling is not 
strictly adiabatic for which energy injection from the inner core 
is needed, or Tbb and Lbb are overestimated by our modelling 
procedure, or underestimated bv lStarUng et alJ (1201 ih . 



Time since the trigger [s] 
10^ 106 



4. Multicolour evolution of SN 201 Obh 

We now present the entire GROND data set, which includes data 
in seven different bands covering the wavelength range from 380 
to 2300 nm. Photometry was corrected for host-galaxy extinc- 
tion computed at the end of Sect. 13.21 In addition, the after- 
glow component derived in Sect. 13.11 was subtracted from the 
data. Assuming a power-law decay for the optical afterglow with 
a = - 1 .3 ± 0.2 computed from the contemporaneous X-ray data, 
a flux and its corresponding uncertainty were derived from the 
power-law model fitted to the early SEDs for each epoch. We 
note however that typically ax + ffopt. The correction for the 
afterglow contribution is more significant at earlier times and in 
redder bands and it has a lower significance at around maximum 
brightness for all bands. 



4. 7. Optical and near-infrared liglit curves 

Figure |5] shows the light curves in the optical g'r'i'z' and the 
NIR JH bandpasses. The g 'r 'i 'z' J light curves show the usual 
pattern of SE SNe: the redder the filter, the later and broader 
the peak. The H band peaks a few days earlier than 7, although 
given the large uncertainties, the time difference is insignificant. 
The band shows no credible detections in any of our observa- 
tions down to limits in the range of 18.4- 18.9 mag (AB system) 
despite our observations covering the expected peak of the SN 
(10.5-25.5 d after the burst). In the following section, we anal- 
yse our data of SN 201 Obh, one of the best-observed GRB-SNe 
to date, and compare the light curve with those of other SNe 
connected to GRBs. 

4.2. Coiour evoiution 

Six GROND colour curves are presented in Fig. |6] to analyse 
the colour evolution of SN 201 Obh. Most colours show drastic 
evolution, where the cooling of the SN photosphere is evident as 
colours become redder. Non-significant variations over time are 
shown in r ' - / ' and J -H colours. 

4.3. Comparison witli SN 1998bw 

To study the luminosity evolution of SN 201 Obh to other SNe, 
we fit SN 1998bw templates d erived fr om the publicly avail- 
able UVBRI photometry (Gala ma et alJ[l998h) to each of our 
GROND filters. The process of compiling these temp lates from 



the o bserved light curve of SN 1998bw is based on IZeh et al 



( 1200 4). Given their modus operandi, the NIR templates are inac- 
curate because they rely on an extrapolation o f the UVBRI data . 
Therefore, the three epochs of JHK^ data in IPatat et al.l (1200 ih 
were used to define the zero points of the flux scale. For all fil- 
ters, it is assume d that the ho st-gala xy extinction of SN 19 98bw 
is equal to zero JPatat et al.ll200ll: Iciocchiatti et al.ll201 ll) . An 
analytical function was then employed to parametrise the SN 




1 10 100 

Time since the trigger [days] 

Fig. 5. Multi-colour light curves corrected for host-galaxy extinction 
and afterglow-subtracted. Filled circles represent detections and arrows 
are upper limits. Solid lines represent the overall fits and dashed lines 
individual components. Fits of SN 1998bw templates are extrapolated 
from K fo + 60 d. For reasons of clarity, light curves were shifted along 
the magnitude axis and the systematical error in Av/_host was not added 
to the error bars. 



1998bw templates. Two of the seven template parameters rep- 
resent the brightness and morphology of the light curve, the fac- 
tors of luminosity k and stretch s, respectively. These were de- 



fined provided that k - 1 and s = 1 for SN 1998bw (IZeh et al 



2004) . In addition, it was necessary to include a delay parameter 
fdeiay that acts by shifting the template linearly in time, i.e., as 
a time offset. The delay parameters that differ from zero do not 
necessarily mean that the onset of the SN is out of phase. 

Given that the early blue emission is inconsistent with the 
templates even after subtraction of the afterglow component de- 
rived in Sect. 13.11 an empirical power-law component was re- 
quired (dashed in Fig. |5]l. The power-law slope was fixed to 
a - -1.3. Leaving a free did not improve the fit and differ- 
ent values of a only negligibly affected the stretch and luminos- 
ity factors derived for the SN component. Table |2] summarises 
the results of fitting the SN 1998bw templates and the empirical 
power-law to the GROND data. The overall fits are shown in Fig. 
|5]for each band using solid lines. 

Luminosity factors listed in the first line of Table |2] which 
include the uncertainty in Ay host, reflect differences between the 
colours of SN 201 Obh and those of SN 1998bw. The g 'band is 
as bright as SN 1998bw and, as a ro ugh compar i son, dis agrees 
with the fainter B-band results from Cano et al. (201 lb) at the 
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Table 2. Fits of SN 1998bw templates to SN 2010bh. 





8' 


r' 


/' 


z' 


J 


H 




SN amplitude (k)^ 


1.10 + 0.12 


0.65 ± 0.05 


0.54 + 0.04 


1.25+0.07 


0.69 + 0.07 


0.63 + 0.08 


< 1.4 


Stretch factor (s) 


0.77 + 0.01 


0.78+0.01 


0.65 + 0.02 


0.73 + 0.02 


0.84 + 0.04 


0.67 + 0.07 




Peak delay (fdeiay)! 


-3.3 


-3.2 


-2.6 


-0.4 


-1.0 








° For afterglow-subtracted and host-extinction-corrected data (A v host = 1.2 + 0.1 mag). 

* Relative to maximum luminosity of SN 1998bw after applying a time offset due to the stretch factor. 



Time since the trigger [s] 
10^ 10** 



-d 
a 

'3 

(0 

s 

P3 




1 10 100 

Time since the trigger [days] 

Fig. 6. Colour curves corrected for host-galaxy extinction and 
afterglow-subtracted. Filled circles represent detections in both filters 
and arrows are upper limits derived from detections in one of the two 
bands. Solid lines represent the colour evolution of SN 1998bw derived 
from the light-curve fits in Sect. 14.31 For clarity, error bars do not in- 
clude the systematic error in Ai/_host and for the last epoch, upper limits 
are not shown. 



5cr confidence level after including the uncertainty in their ex- 
tinction determination. Our r 'i 'luminosity factors are also larger 
than previously reported from opti cal photometry of SN 2010bh 
(k » 0.4-0.5: ICano et al.ll201 Ibl) mainly owing to the use of 
a different extinction correction. There is certainly some intrin- 
sic blue excess in the g ' band, although the difference from the 
r'band of 40% has a rather large uncertainty of 13%. Higher 
luminosities at bluer wavelengths indicate that the photosphere 
has a higher temperature at its peak luminosity than that of SN 
1998bw. Given that the expanding atmospheres of SNe cool 
down over time and SN 2010bh peaks about a week earlier than 
SN 1998bw in g ', it is a natural conclusion that it must be hotter 



at maximum brightness. It is also possible that a g -band peak 
as early as 8 d after the trigger still contains a non-negligible 
contribution from the cooling shock breakout. We also note that 
the z -band luminosity factor is much larger than in the other 
bands; it is brighter than SN 1998bw at the 3cr confidence level. 
Af ter inspection of the spectra around maximum light presented 
bv lBufano et al.l (1201 ih . the Ca ii /18579 emission Une was found 
to contribute roughly 2 - 4% to the continuum flux integrated in 
the z'-band sensitivity range (8254-9528 A). The contribution 
of this spectral feature is not enough to account for the excess 
of 45% in the z 'band relative to the infrared. The discrepancy is 
therefore attributed to the systematical u ncertainties introduced 
when extrapolating the UBVRI data from Galama et al. ( 1998bl) 
to construct the templates of SN 1998bw. In summary, the lumi- 
nosity of SN 2010bh is a factor 0.5 - 0.7 fainter than SN 1998bw 
in optical r'/' bands and 0.6-0.7 in NIR JH bands. We also 
note that at late times SN 2010bh fades more rapidly in the op- 
tical than SN 1998bw did, although the fluxes if the templates 
were extrapolated after ^ fo + 60 d. 

The stretch factors listed in the second line of Table |2]range 
from 0.6 to 0.8 and are at the low end of the GRB-SNe distribu- 
tion (see figure 5 of Ferrero et al. 20061). T hese are also consis- 
tent with the findings bv lCano et al.l(l201 lb!) for SN 2010bh. The 
optical light curves are wider than the stretch factor suggests, or 
in other words, the stretch factor predicts later peak times for the 
optical light curves. To account for this, fdeiay shifts the templates 
to earlier times by about three days in g 'r 'i ', which means that 
the optical light curves of SN 2010bh peak even earUer than the 
stretch factor suggests. After the inclusion of the delay parame- 
ter fdeiay, the strctch factors now solely represent the width of the 
light curves and not the time of the maximum brightness com- 
pared to SN 1998bw. Moreover, peak times of roughly 8 - 9 d in 
the g 'r 'bands (equivalent to the V in the Johnson filter system) 
corresponds to the earliest and fastest light curve s of GRB-SNe 
observed to date (see table 2 of Richardson 12009V 

Offsets in peak time, in our case represented by fdeiay, can 
be explained as the result of delayed black-hole formation 
(IVietri&Stellalll999h . n this scenario, the XRF might be trig- 
gered by the core collapse of the progenitor to a neutron star, 
soon after which accretion holds. The supernova would then oc- 
cur after the further collapse of the neutron star into a black hole. 
The delay could b e of the order of months or years or perhaps as 
short as hours (see lZeh et al ]|2004 . Nevertheless, the time delay 
here is negative and it is much more plausible that the compar- 
ison of light-curve morphology to SN 1998bw may need more 
than two parameters to be accurate. 

Fits of SN 1998bw templates were also employed to study 
the colour evolution of SN 2010bh in detail as presented in Fig. 
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|6l where the curves are shown without the empirical afterglow 
component. From the first two data points, it is clear from the 
g'-r' and g'-z' colours that there is a blue component that can- 
not be modelled by the templates and is interpreted as the shock 
breakout in Sect. 13.11 From these two colours, it is also possi- 
ble to see that SN 2010bh becomes red faster than SN 1998bw 
did. The r'-i' and J -H colours remain roughly constant, which 
shows that the changes occur on a broader wavelength scale. The 
standard colour evolution from blue to red is shown by r' - z' 
and r' - J, which at late times evolve bluer and redder than SN 
1998bw templates, respectively. 

4.4. Comparison with SN 2006aj 



We no w compare results from optical data by iFerrero et al 
who used the same technique and templates to get the 
luminosity and stretch factors for SN 2006aj. They computed 
luminosity factors in the range 0.62-0.76, which are approxi- 
mately in the same range as SN 2010bh without considering the 
measurements in g 'z '■ Their stretch factors ranged from about 
0.62 to 0.69, which makes the SN 2010bh optical light curves 
wider than those of SN 2006aj. In contrast to the definition of 
the stretch factor, namely that earlier peak times tend to corre- 
spond to a narrower light curve, peak times are earlier in the 
case of SN 2010bh. Whilst the BV photometry for SN 2006aj 
peaked roughly 9 and lid after the burst, respectively, the g'r' 
photometry for SN 2010bh peaks approximately 8 and 9 d after 
trigger, respectively. This supports the statement that SN 2010bh 
has evolved more rapidly than any other GRB-SNe, given that 
SN 2006aj had the fastest evolution until th e discovery of SN 
2010bh. The JH ligh t curves of SN 2006aj dCobb et al.ll2006l 
Kocevski et al.l 120071) are scaled to the luminosity distance of 
SN 2010bh for comparison. After the host-extinction correction, 
SN 2010bh turns out to be as bright as SN 2006aj in the NIR 
as well. Hence, SN 2010bh is similar to SN 2006aj in terms of 
light-curve shape and luminosity. 



5. Bolometric light curve 

The bolometric light curves of SNe are an essential tool for ex- 
amining global luminosity features and enable us to compare 
with other SNe and theoretical models. However, it is difficult 
to obtain such a light curve because of the limited information 
at UV, infrared, and radio wavelengths. Only a ^Mai/-bolometric 
light curve can be constructed by using a broad spectral coverage 
to derive a total flux that is then used as a proxy of the bolomet- 
ric flux. To accomplish this task, we employed the wavelength 
range covered by our g'r'i'z' JH filters, i.e., from 350 to 1800 
nm. 

In the case of SN 2010bh, there are no UV constraints. In the 
following description, no attempts at correcting for the UV flux 
were made. By using the afterglow-subtracted g'r'i 'z' JH pho- 
tometry corrected for host-galaxy extinction, monochromatic 
fluxes for each bandpass were derived. Sets of three bandpasses 
were defined to interpolate their corresponding monochromatic 
fluxes using the Simpson's rule. The second-degree polynomial 
result of the interpolation was then integrated over frequency in 
the range of each set of bandpasses. Finally, the total flux in 



the range from 350 to 1800 nm was determined by adding up 
the integrated fluxes of all sets of bandpasses. The total flux is 
transformed to a quasi-bolometric luminosity using a distance 
of 240 ± 17 Mpc to the host galaxy of SN 2010bh (see Sect.[3j] 
for more details). No attempts of extrapolation beyond the lim- 
its of the g ' and the H bandpasses were made. Corrections for 
the NIR flux at late times were found to be the most significant. 
The data at fo + 30.7 d (rest frame) were corrected for the //-band 
non-detection by assuming that the fraction of //-band flux com- 
pared to the total bolometric flux remains constant at 8% starting 
from fo -H 24. 1 d. Similar corrections were performed for JH non- 
detections at 38.6 and 78.8 d after the burst under the assumption 
of a constant JH flux fraction of 27% at fo -i- 30.7 d. Results are 
shown in Fig. |7] along with quasi-bolometric hght curves from 
SB and other GRB-SNe. 

The analysis of the quasi-bolometric light-curve morphology 
yields a peak luminosity of 4.3x10^^ erg s"' at about 8 d after 
the trigger (equivalent to Mi,oi ~ -17.87), i.e., approximately 
two times fainter and six days sooner than for S N 1998bw. 
Our lum inosity is 16% higher than that computed bv lCano et al. 
(1201 Ibl) . although consistent to within our 1 1 % of uncertainty 



Whilst the early peak of SN 2010bh correlates with its narrow 
ness and low luminosity, this is the case for neither the entire 



GRB-SNe sample nor the local sample of SB SNe dZeh et al 



2004; Richardson et al. 2006; Richard sonir2009l) . Th e morphol- 



ogy of the light curve is similar to that of SN 2006aj (IPian et al 
2006), although 21% faint er. The peak time also resembles that 



of the type-Ic SN 19941 ([Richmond et al.lll996l) . although SN 
2010bh has a much wider light curve, which is 77% brighter at 
maximum. In terms of peak luminosit y, SN 2010bh is similar to 
the broad-lined Ic SN 2009bb (Pignat a etalll201lh . It also un- 
derwent the most dramatic late-time decay in the sample, which 
implies that its envelope became rapidly optically thin to y-rays. 
The last statement is supported by the extremely high expansion 
velo cities measured for SN 2010bh of the order of 30,000 km 
s" dChornock et al. I l2010ah . Another clear feature is the sudden 
decrease in luminosity at around fo + 30 d, which contrasts with 
the smooth decay in the comparison SNe at similar stages. This 
indicates either that the atmosphere becomes rapidly optically 
thin or that the assumption of a constant NIR contribution after 
w fo H- 3 1 d underestimates the flux in the JH bands. 

5. 1 . Physical parameters of tlie explosion 



We followed the approach described in IValenti et al.l (l2008h 
to derive the physical quantities that characterise the explo- 
sion, i.e., we modelled the early and late light curves sepa- 
rately. The early-time phase corresponds to the photo spheric 
regime for which the analytical model developed by lArnett 
( 119821) has been adopted, initially used for SNe la and adapted 
to SB SNe (e.g..lTaubenbergeretal.ll2006t IValenti etalj|2008 : 
IPignata et ai]|201 ll: iBenetti et al.ll201 ih . At late stages, the at- 
mosphere becomes nebular, i.e., optically thin, and the emitted 
luminosity is powered by the energy deposition of: (1) y-rays 
from ^^Co decay, (2) y-rays from electron-positron annihilation, 
and (3) the kinetic energy of the positron s (see appendix A in 
Valenti et a02008h . However. iMaeda et all (l2003l) noted that the 



two-component configuration leads to inconsistencies between 
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Fig. 7. Quasi-bolometric light curve of 
SN 2010bh produced by using GROND 
g'r'i'z' JH filters (black open circles) 
in the rest frame. For clarity, the uncer- 
tainties in the host-galaxy extinction 
are not included in the error bars. A 
single black continuous line represents 
the best-fit model (see Sect. 15. U . Early 
and late components of the model 
are smoothly joined at Iq + 30 d. SE 
and other GRB-SNe quasi-bolometric 
light curves have been plotted as a 
comparison sample: SN 20031w (ORB 
03 1 203 ; Malesani et al. 20 04J, SN 
1998bw (ORB 980425; (Galarna et all 
il998bl). SN 2 006ai (ORB 060218; 
iPian et alj|2006r). the broad -lined Ic SN 
2009bb ("Pignata et al.'201l(), the ty pe-Ic 
SN 19941 (Richmond et all 1 1996b . the 
type Ibc SN 2008D (XRO 080109; 
Modjaz et al] 120091 ; ISoderberg et"ai] 
20081), and the type-Ic SN 2002ap 



2UU«I). and ttie tvpe-lc ^N 2UU2ap 
dOal-Yam et al.ll2002l ; lFoievetai]|2003l ; 



lYoshii et alJl2003l)~ 



the parameters derived from fitting the early and late light curves 
of SNe Ic, caused mainly by the model limitations in varying the 
y-ray trapping over time. To enable low and high y-ray trapping 
at early and late times, respectively. iMaeda et al.l (l2003b divided 
the ejecta into a high-density inner region and a low-density 
outer region. The emission from the outer region dominates the 
total emission in the optically thick regime at early times, and 
that from the inner region, which has a higher y-ray opacity, 
dominates in the nebular phase at late times. Here, we use the 
same procedure to model the g'r'i 'z 'JH quasi-bolometric light 
curve of SN 2010bh. 

Given the model explained above, a total of four free pa- 
rameters were used to fit the quasi-bolometric light curve of SN 
2010bh: the total mass of ^^Ni produced in the envelope Mni, the 
total ejecta mass Mgj, the fraction of mass in the inner component 
/m, and the fraction of kinetic energy in the inner component /£. 
The kinetic-energy-to-ejected-mass ratio of the outer region was 
fixed by usi ng its correla tion with photospheric velocity at peak 
luminosity (lArnettlll982h 



3 2£k.< 



5 M 



(3) 



ej.out 



This expression assumes that the density of the ejecta is homo- 
geneous and that the inner component does not contribute to 
the emitted luminosity in the optically-thick regime. Since the 
photospheric velocity was not available directly from observa- 
tions, the velocity measured by fitting P-Cygni line profiles was 
used as a proxy of Uph. However, the envelope layer where the 
blue-shifted absorption line forms does not necessarily coincide 
with the position of the photosphere, as found when measuring 
different expansion velocities from absorption lines of different 
species. The spre ad can amount to several hundreds km s ' (see 
Jones et aI.lEo09l for an example of type-II SNe). For the spec- 



tra of SN 2010bh, IChornock et all feOIOal) obtained velocities 
of about 35,000 and 26,000 km s"' from the Si ii .16355 feature 
roughly 21 and 6 d after the burst, respectively. Since there is 
no measurement at the time of maximum light, which is about 
8 - 9 d after the trigger, a range of expansion velocities was used 
in the modelling. Assuming that the photosphere lies at deeper 
layers than those where lines are formed and recedes in mass ex- 
posing deeper and slower layers, in the modelling we employed 
photospheric expansion velocities of 2.5, 2.8, and 3.1 xIO'* km 
s Other physical and mathematical quantities such as opac- 
ity and integration constants were chosen to be the same as in 
Cano etal.l(20IIbl). 



Our fitting procedure consisted of two steps. First of all, we 
modelled the data around maximum luminosity {5 < t - to < 
30 d) using the Arnett's model, assuming that only the outer 
component contributes to the total luminosity at this stage. We 
then obtained MNi.out and Mej,out- Secondly, we modelled the late- 
time data (t - tg > 30) using the nebular-phase components of 
Valenti et al.l (120081) assuming that both the inner and outer re- 
gions contribute to the total emitted luminosity. Here, we fixed 
A/Ni,oiit and Mej.out to the values obtained in the first step of the 
fitting procedure and only /m and were allowed to vary. 

From the above modelling scheme, our best-fit parameters 
were MNi.out = (0.135 + 0.001)Mo, Mej.out = 2.37-2.90 Mq, 
/m = 0.36 + 0.04, and = 0.1 1 -0.15 for the three selected ex- 
pansion velocities at the photosphere. All of these results com- 
bined together provided the total masses and energy of the ex- 
plosion listed in Table [3] (statistical eiTors only). The ^^Ni mass 
is independent of the chosen expansion velocity at the photo- 
sphere, Mni - (0.21 + 0.03)Mq, given that it is proportional 
to the luminosity. In contrast, since the ejected mass changes 
significantly, the weighted mean and the RMS of the results of 
the three models in Table |3] were employed to compute a fi- 
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Table 3. Physical parameters of the explosion with a varying expansion velocity. 



Vph 

[kms"'] 


A^Ni.out 

[Mo] 


-^ej,out 

[Me] 


Im 


fE 


Mm 
[Me] 


[Me] 


[10^2 erg] 


25,000 
28,000 
31,000 


0.135 + 0.001 
0.135 + 0.001 
0.135 + 0.001 


1.52 + 0.05 
1.70 + 0.05 
1.87 + 0.06 


0.36 + 0.02 
0.36 + 0.04 
0.36 + 0.02 


0.15 + 0.03 
0.13 + 0.07 
0.11+0.02 


0.21+0.02 
0.21+0.03 
0.21+0.02 


2.37+0.10 
2.64 + 0.14 
2.90 + 0.12 


1.85+0.09 
2.52 ± 0.22 
3.37+0.13 



nal value of Mej = (2.60 ± O.23)M0. The total kinetic energy 
and energy fraction were derived in the same way, implying that 
£k = (2.4 + 0.7) X 10^2 erg and /e = 0.12 ± 0.02. 

Using comparable independent data, Cano et al. ( 201 Ibl) fol- 
lowed a similar procedure to derive physical quantities from 
quasi-bolom etric data. W hilst our values for are consistent 
with those in lCano et al. fcol Ibl) to within 1.5cr, they found that 
Mm = (0.10 ± O.O1)M0 for SN 2010bh, which is two times 
lower than our value. The discrepancy affects the determination 
of the kinetic energy as well, which is connected in Eq.[3]to Wph 
and Mej. The causes of these inconsistencies are: (1) the host- 
galaxy extinction employed, which is E (B - V )host = 0.39 ± 0.03 
mag in our case and 0.18 ± 0.08 mag in lCano et al. ( 201 Ibl) : (2) 
the inclusion of an inner component hidden at optically thick 
stages, which increases the ^^Ni mass by a factor of /m', and (3) 
the choice of different expansion velocities of the photosphere. 
Furthermore, even when including the uncertainties in the host- 
ext inction determinatio n, which are cthj^^ ~ 0.02 both here and 
for Cano et al. ( 201 Ibl) . the discrepancy persisted. Nevertheless, 
evidence for the existence of a dense inner layer of 0.94+0. 15M0 
is provided and supported by the sub-luminous post-maximum 
phase of SN 2010bh, which is indicative of a high trapping of 
y-rays. 

The large amount of ""^Ni produced is consistent with the 
class of SN associated with GRBs, but also matches the value de- 



rived for the highly energetic type-Ic SN 2004aw (IValenti et al 



12008.) . which in contrast had a much broader light curve that 
varied more slowly with time. Most interestingly, SN 2010bh re- 
sembles SN 2006aj in terms of light-curve shape. The ^^Ni mass 



produ ced in the explosion is practically the same (IMazzali et al 



20061) . whilst the ejected mass is 20% higher for SN 2010bh and 
the total kinetic energy is significantly different (~ 2 x 10^' erg 
in the case of SN 2006aj). This ensures that SN 2010bh is re- 
markable in terms of expan sion velocity, which ha s been one of 
the greatest ever measured (IChornock et al ]l2010al) . 

We note that the end of the optically-thick regime and the 
beginning of the nebular phase cannot be accurately defined. 
The rule of thumb is that the photospheric phase ends 30 d after 
explosion at the earliest and that the nebular phase starts 60 d 
after the explosion at the latest. Since Arnett's model fits our 
data relatively well until day 30, it is defined as the end of the 
optically-thick era. Nevertheless, we are unable to clearly es- 
tablish whether the period between 30 and 60 d after the burst 
corresponds to the nebular phase already, as the ejecta are most 
probably neither completely thick nor sufficiently thin. 

The model reproduces the luminosity at maximum bright- 
ness and thereafter, despite underpredicting the luminosity dur- 
ing the rising phase. This is because the light curve is domi- 
nated by shock-breakout emission at early stages. The same is 



obtained in Sect. |4| when trying to fit templates of SN 1998bw 
to our multicolour light curve, where an extra component for 
the shock breakout is required. The light curves of SN 2010bh 
peak so early that an even more rapid evolution than observed 
is expected after maximum luminosity. This early-wide peak di- 
chotomy could be explained by the ^^Ni distribution in the enve- 
lope. If there were more ^^Ni produced near the surface, the peak 
would be early and wide, i.e., the SN would rise more slowly but 
much earlier In contrast, if the ^^Ni were concentrated towards 
the centre, it would lead to a much later peak, although sharper 
rise, producing a narrower sh ape of the maximum (M. Bersten, 
2010, private communication : lNomoto et al hold) . 



6. Conclusions 

Spanning a time range from 12 hours to 83 days after the trig- 
ger and covering from 190 to 2300 nm in wavelength (see 
Sect. IS, we have presented UV/optical/NIR photometric data 
of XRF 100316D/SN 2010bh. Given the results inti'oduced and 
discussed in Sect. [3] ID andjSj we have drawn the following con- 
clusions: 

- Broad-band SEDs at early times demonstrate the existence of 
red and blue components identified as synchrotron emission 
from the XRF afterglow and the cooling envelope after shock 
breakout, respectively. 

- A significant amount of dust along the line of sight through 
the host galaxy (Ay host = 1.2 + 0.1 mag) is consistent with 
the two-component SED model and agrees with the faint de- 
tections at UV wavelengths. 



- By co mparing with earlier X-ray results from lStarling et al 



(1201 ll) . we have demonstrated that the temperature of the 
blackbody component decreases with time, which is consis- 
tent with a scenario of a cooling expanding atmosphere. 

- By performing an additional analysis of the thermal compo- 
nent and the earlier X-ray measurements we have measured 
expansion velocities that are consistent with SN expansion 
and an initial apparent emission radius of 7 x 10^' cm. This 
radius is slightly larger than the size of WR stars, which are 
the most likely GRB-SNe progenitors. If a WR star were 
the progenitor of XRF 1003 16D/SN 2010bh, then the initial 
radius could indicate that there was a massive dense stellar 
wind surrounding the progenitor 

- Our multicolour light curves after host-galaxy correction and 
subtraction of the afterglow component, have peak r'i' lu- 
minosities of about 0.5-0.7 times that of SN 1998bw and 
consistent with those of SN 2006aj. Similarly the NIR lumi- 
nosity at maximum is as bright as SN 2006aj and 0.6-0.7 
times that of SN 1998bw. The excess in the ^'band indi- 
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cates that SN 2010bh has a hotter photosphere than that of 
SN 1998bw at the time of maximum brightness. 

- We have found that SN 2010bh is the most rapidly evolv- 
ing GRB-SNe to date, reaching maximum luminosity 8-9 
days after the burst in the bands. At late times, it also 
fades more rapidly than SN 1998bw showing redder colours 
as well. This behaviour is also evident in the bolometric light 
curve, which decays faster than for any SN in the comparison 
sample. 

- The physical parameters of the explosion are derived by 
means of the quasi-bolometric light curve constructed from 
our g'r'i'z 'JH pho tometry. The modelling is performed us- 
ing Arnett's model (lArnettill982r) for data around peak and 
standard y-ray deposition at later times. A high-density in- 
ner component with roughly 26% of the total mass is re- 
quired to reproduce the flux ratio between maximum lumi- 
nosity and tail. The total mass of ^^Ni produced in the en- 
velope is Mm = 0.21 + O.O3M0, which precisely matches 
the value derived for SN 2006aj, whilst the total ejecta mass 
of Mej = 2.6 + O.2M0 exceeds the value for SN 2006aj by 
20%. However, the kinetic energy turns out to be higher at 
Ek = (2.4 + 0.7) X 10^2 erg, making SN 2010bh the second 
most energetic GRB-SN after SN1998bw. 

The association between XRF 100316D and SN 2010bh is 
particularly interesting, since for the second time the cooling of 
the shock breakout has been detected in a GRB-SN. It is also 
unique in revealing a hot component that possibly contributes 
even at blue maximum brightness, one of the largest host-galaxy 
extinctions measured for this kind of transient, and the fastest 
rise among GRB-connected SNe. 
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Appendix A: Optical and near-infrared photometry 



Table A.l. GROND photometry of field stars used for relative photometry. 



R.A.(J2000) Dec.(J2000) g' 


r' 


i' 


z' 


J 


H 



















07 
07 

07 
07 
07 
07 



10:30.10 
10:34.37 

10:28.77 
10:31.17 
10:29.97 
10:30.98 



-56:14:58.1 
-56:15:58.1 
-56:14:30.9 
-56:14:45.8 
-56:16:24.4 
-56:14:26.3 



17.981(16) 
18.092(16) 

16.221(13) 
16.232(12) 
17.703(13) 
16.804(13) 



17.423(10) 
17.365(10) 

15.831(10) 
15.535(09) 
16.435(10) 
16.391(10) 



17.246(12) 
17.158(11) 

15.742(11) 
15.305(10) 
15.902(10) 
16.267(11) 



17.065(16) 
16.973(16) 

15.596(15) 
15.094(15) 
15.570(15) 
16.117(16) 



16.587(12) 
16.477(12) 

15.224(07) 
14.606(06) 
14.875(06) 
15.733(09) 



16.737(13) 
16.646(13) 

15.404(08) 
14.642(07) 
14.864(07) 
15.949(09) 



16.964(31) 

15.933(15) 
15.083(11) 
15.216(11) 
16.374(22) 



Notes. All magnitudes are in the AB system and corrected for neither Galactic nor host-galaxy extinction. Errors are statistical only. 



Table A.2. GROND photometry of XRF 1003 16D/SN 2010bh after image subtraction. 



Time Interval 
[d] after the trigger 




g' 


r' 




i' 




z' 


J 




H 


0.49023- 


0.49470 


20 


75(04) 


20.90(04) 


20 


87(05) 


20 


77(07) 


20.69(14) 


> 


20.04 


0.56046- 


0.60240 


20 


81(03) 


20.91(03) 


20 


94(04) 


20 


76(04) 


20.71(26) 


> 


20.67 


1.47634- 


1.49123 


20 


75(04) 


20.90(04) 


20 


59(05) 


20 


51(05) 


20.47(22) 


20 


69(33) 


1.69356- 


1.70396 


20 


92(04) 


20.77(03) 


20 


69(05) 


20 


65(06) 


20.40(32) 


> 


20.09 


2.49804- 


2.51041 


20 


70(04) 


20.58(04) 


20 


50(05) 


20 


20(04) 


19.99(17) 


> 


20.09 


3.50557- 


3.52164 


20 


67(03) 


20.40(03) 


20 


24(03) 


19 


94(03) 


19.72(17) 


19 


94(21) 


4.48027- 


4.48963 


20 


43(04) 


20.21(03) 


20 


15(04) 


19 


78(04) 


20.00(19) 


20 


22(23) 


5.46324- 


5.47229 


20 


33(04) 


20.11(04) 


20 


04(04) 


19 


61(04) 


19.93(19) 


19 


92(23) 


7.47679- 


7.48174 


20 


26(06) 


19.96(04) 


19 


87(06) 


19 


58(06) 


19.69(19) 


20 


20(28) 


10.47338- 


10.48259 


20 


38(03) 


19.91(03) 


19 


98(04) 


19 


37(03) 


19.53(14) 


19 


85(25) 


11.47550- 


11.47994 


20 


50(03) 


20.01(03) 


19 


99(04) 


19 


33(03) 


19.49(14) 


19 


94(18) 


14.47912- 


14.49638 


20 


71(04) 


20.07(03) 


20 


08(03) 


19 


38(03) 


19.51(12) 


19 


67(22) 


15.48101- 


15.48997 


20 


72(05) 


20.07(04) 


20 


16(05) 


19 


32(04) 


19.59(13) 


19 


85(19) 


19.52885- 


19.53123 


21 


33(10) 


20.36(05) 


20 


39(07) 


19 


51(05) 


19.70(19) 


> 


19.71 


21.48630- 


21.49543 


21 


52(06) 


20.48(04) 


20 


48(04) 


19 


65(03) 


19.64(13) 


19 


95(19) 


25.49771- 


25.51874 


21 


91(05) 


20.78(03) 


20 


66(03) 


19 


87(03) 


19.67(13) 


20 


30(21) 


32.48587- 


32.50697 


22 


78(10) 


21.42(04) 


21 


28(05) 


20 


39(04) 


20.07(19) 


> 


20.41 


40.46892- 


40.50042 


> 


23.11 


21.83(06) 


21 


74(07) 


20 


81(05) 


> 20.34 


> 


20.43 


83.43058- 


83.45161 


> 


23.95 


23.66(23) 


> 


22.91 


22 


78(25) 


> 20.95 


> 


20.90 



Notes. The template image for subtraction was taken 234 days after burst. The time reference is to = 55271.531 13 MJD, the date of the GRB 
trigger. All magnitudes are in the AB system and corrected for neither Galactic nor host-galaxy extinction. Upper limits are all at the 3a 
confidence level. Errors include the systematics of the photometric calibration. 



